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There is nowadays considerable interest on ways to quantify the dynamics of research activities, 
in part due to recent changes in research and development (R&D) funding (lH9J. Here, we seek to 
quantify and anal yze university research activities, and compare their growth dynamics with those 
of business firms |10h14J . Specifically, we analyze five distinct databases, the largest of which is a 
National Science Foundation database of the R&D expenditures for science and engineering of 719 
United States (US) universities for the 17-year period 1979-1995. We find that the distribution of 
growth rates displays a "universal" form that does not depend on the size of the university or on the 
measure of size used, and that the width of this distribution decays with size as a power law. Our 
findings are quantitatively similar to those independently uncovered for business firms (RMlji, and 
consistent with the hypothesis that the growth dynamics of complex organizations may be governed 
by universal mechanisms. 
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In the study of physical systems, the scaling properties 
of fluctuations in the output of a system often yield in- 
formation regarding the underlying processes responsible 
for the observed macroscopic behaviour JlaJlq]. Here, 
we analyze the fluctuations in the growth rates of uni- 
versity research activities, using five different measures 
of research activity. The first measure of the size of a 
university's research activities that we consider is R&D 
expenditures. The rationale for using R&D expenditures 
as a measure of research activity is that research is an 
expensive activity that the university finances with ex- 
ternal support. 

We first analyze a database containing the annual R&D 
expenditures for science and engineering of 719 US uni- 
versities for the 17-year period 1979-1995 (« 12,000 
data points). The expenditures are broken down by 
school and department. The annual growth rate of R&D 
expenditures is, by definition, g(t) = log[S(t -I- 1)/S(t)], 
where S(t) and S(t + 1) are the R&D expenditures of 
a given university in the years t and t + 1 respectively. 
We expect that the statistical properties of the growth 
rate g depend on S, since it is natural that the fluctua- 
tions in g will decrease with S. Therefore, we partition 
the universities into groups according to the size of their 
R&D expenditures (Fig. ||a). Figure |]b suggests that 
the conditional probability density, p(g\S), has the same 
functional form, with different widths, for all S. 

We next calculate the width a(S) of the distribution 
of growth rates as a function of S. Figure [lp shows that 
a(S) scales as a power law 



a(S) ~ S-P, 



(1) 



with [3 = 0.25 ± 0.05. In Fig. |l|d, we collapse the scaled 
conditional probability distributions onto a single curve. 
To test if these results for the dynamics of R&D ex- 
penditures are valid for other measures of research ac- 
tivity, we next analyze another measure of a university's 



research activities, the number of papers published each 
year |18|-f20|]. We analyze data for the 17- year period 
1981-1997 from the US University Science Indicators 
pl[ , which records the number of papers published by 
the top 112 US universities (~ 1,900 data points). We 
find that the analog of Fig. fy holds. Particularly strik- 
ing is the fact that the same exponent value, j3 = 1/4, 
is found (Fig. 0a) and that the same functional form of 
p(g\S) is displayed (Fig. |b). 

Next, we consider as a measure of size the number of 
patents issued to a university j p2[ . We "manually" re- 
trieve from the webpages of the US Patent and Trade- 
mark Office's database Q the number of patents is- 
sued to each of 106 universities each year of the 22-year 
period 1976-1997 (« 2 300 data points). We confirm 
that the analog of Fig. [j] holds, with the same exponent 
value, f3 — 1/4 (Fig. Ha), and the same functional form 
ofp(ff|S),Fig.g>. 

To test if our findings hold for different academic sys- 
tems, we analyze two databases on research funding of 
English p3] and Canadian p5| universities. The same 
quantitative behavior is found for the distribution of 
growth rates and for the scaling of a, with the same ex- 
ponent value (Fig. ||a) and the same functional form of 
p{g\S), Fig. 0b. Thus, the analysis of all five databases 
confirms that the same quantitative results hold across 
different measures of research activity and academic sys- 
tems. 

We next address the question of how to interpret our 
empirical results. We start with the observation that re- 
search is an expensive activity, and that the university 
must "offer" its research to external sources such as gov- 
ernmental agencies and business firms. Thus, an increase 
in R&D expenditures at university A and a decrease at 
university B implies that the funders of research increas- 
ingly choose their research from university A as opposed 
to university B [jy . This qualitative picture parallels the 



competition among different business firms, so it is nat- 
ural to enquire if there is quantitative support for this 
analogy between university research and business activi- 
ties. To quantitively test this analogy, we note that the 
results of Fig. [j| ar e remarkably similar to the results 
found for firms |l^Jl4[| and countries |2q-|2q|. We plot 
in Fig. ||c the scaled conditional probabilities p(g\S) for 
countries, firms and universities, and find that the distri- 
butions for the different organizations fall onto a single 
curve. 

There is, however, one difference: For firms and coun- 
tries, we find (3 « 1/6, while for universities, (3 « 1/4. 
We can understand this difference using a model for or- 
ganization growth p9fl . In the model, each organization 
— university, firm, or country — is made up of units. The 
model assumes these units grow through an indepen- 
dent, Gaussian-distributed, random multiplicative pro- 
cess with variance W 2 . Units are absorbed when they 
become smaller than a "minimum size" , which is a func- 
tion of the activity they perform. Units can also give 
rise to new units if they grow by more than the mini- 
mum size for a new unit to form. The model predicts 
f3 = W/[2(W + £>)], where V is the width of the distri- 
bution of minimum sizes for the units p9[ . For firms, the 
range of typical sizes is very broad — from small software 
and accounting firms to large oil and automobile firms — 
suggesting a large value of T>. On the other hand, for 
universities, the range of typical sizes is much narrower, 
suggesting a small value of T> and implying a larger value 
of f3 than for business firms. This is indeed what we 
observe empirically. 

Business firms are comprised of divisions and universi- 
ties are made up of schools or colleges, so it is natural to 
consider the internal structure of these complex organi- 
zations pOfl . We next quantify how the internal structure 
of a university depends on its size by calculating the con- 
ditional probability density p(£\S) to find a school of size 
£ in a university of size S (Fig. Ha). The model predicts 
that p(C\S) obeys the scaling form pg] 



p(t\S)~S- a f(Z/S a ), 



(2) 



where f(u) ~ u T for u <C 1, and f(u) decays as a 
stretched exponential for «»1. We find r = 0.37 ±0.10 
(Fig. |b), and a = 0.75 ± 0.05 (Fig. |c). We test the 
scaling hypothesis (0) by plotting the scaled variables 
p(£\S)/S~ a versus £/S a . Figure ||b shows that all curves 
collapse onto a single curve, which is the scaling function 

/(«)■ 

Equation (0) implies that the typical number of schools 
with research activities in a university of size S scales as 
S 1 ~ a , while the typical size of these schools scales as S a . 
Hence, we can calculate how a depends on S, 



a{S) ~ iS 1 "*)- 1 ' 2 W(£). 



(3) 



7 = 0.16 ± 0.05. Substituting into (||) and remember- 
ing that the typical size of the schools is S a , we obtain 
a(S) ~ (S 1 - a )- 1 / 2 (S a )-" i , which leads to the testable 
exponent relation 



P = 



1-a 



- OTj . 



(4) 



For a sw 3/4 and 7 w 1/6, Eq. (g) predicts (3 « 1/4, 
in surprising agreement with our empirical estimate of (3 
from the five distinct databases analyzed (Fig. 0a). 

Our results are consistent with the possibility that 
the statistical properties of university research activities 
are surprisingly similar for different measures of research 
activity and for distinct academic systems. Moreover, 
our findings for university research resemble those inde- 
pendently found for business firms |10|-|l4| and countries 
p6| p8[ . One possible explanation is that peer review, 
together with government oversight, may lead to an out- 
come similar to that induced by market forces, where the 
analog of peer-review quality control may be consumer 
evaluation, and the analog of government oversight may 
be product regulation. 



In order to determine a, we first find the dependence 
of W on £. Figure |d shows that W ~ £~ 7 with 
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FIG. 1. Growth dynamics of research activities at uni- 
versities, a) Histogram of the logarithm of the annual R&D 
expenditures of 719 US universities for the 17-year period 
1979-1995, expressed in 1992 US dollars. Here, 5" denotes the 
R&D expenditures detrended by inflation so that values for 
different years are comparable. The bins were chosen equally 
spaced on a logarithmic scale with bin size 0.5. The line is 
a Gaussian fit to the data, which is a prediction of Gibrat's 
theory [ |10|Jl2| . b) Conditional probability density function 
p(g\S) of the annual growth rates g. For this plot the entire 
database is divided into three groups (depicted in (a) by dif- 
ferent shades), c) Standard deviation cr(S) of the distribution 
of annual growth rates of as a function of S. The straight line 
is a power law fit to the the data, and its slope gives the expo- 
nent /3 — 0.25 ± 0.05. d) Scaled probability density function 
p(g\S)/a~ 1 (S) plotted against the scaled annual growth rate 
(g — g)/a(S) for the three groups defined in (b). Note that 
the scaled data collapse onto a single curve. 
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FIG. 2. Robustness of empirical findings for the distri- 
bution of growth rates, a) Standard deviation a(S) of the 
distribution of annual growth rates for different measures of 
research activities and different academic systems from the 
data in the five distinct databases analyzed: (i) the number 
of papers published each year at 112 US universities, (ii) the 
number of patents issued each year to 106 US universities, (iii) 
the R&D expenditures in US dollars of 719 US universities, 
(iv) the total amount in Canadian dollars of the grants to 60 
Canadian universities, and (v) the external incomes in British 
pounds of 90 English universities. It is apparent that for all 
measures and all academic systems analyzed, we find a power 
law dependence — with the same exponent j3 ~ 1/4. The 
values of a for the different measures were shifted vertically 
for better comparison of the estimates of the exponents, b) 
The distribution of annual growth rates, scaled as in Fig. Fyd, 
for the five databases. We show the distribution of growth 
rates for 2 different groups, obtained in a way similar to that 
described in Fig. |jb, for each of the five measures. The data 
appear to collapse onto a single curve, suggesting that the 
different measures have similar statistical properties, c) The 
distribution of scaled annual growth rates for different organi- 
zations: R&D expenditures of US universities, sales of firms, 
and GDP of countries. The data collapse onto a single curve 
suggesting that the scaled distributions have the same func- 
tional form. 
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FIG. 3. Statistical analysis of the units forming the in- 
ternal structure of a university, the schools, a) Conditional 
probability function p(£\S) of finding a school of size £ in 
a university of size S. To improve the statistics, we parti- 
tion the universities by size into four groups, b) To illustrate 
the scaling relation (bl), we plot the scaled probability den- 
sity p(£\S)/S~ a versus the scaled size of the school £,/S a . In 
agreement with (Q), we find that the scaled data fall onto a 
single curve, c) Scaling of the typical size of a school in a 
university of a given size for different university sizes. The 
data obey a power law with exponent a = 0.75 ± 0.05. d) 
Standard deviation W of the distribution of growth rates of 
schools versus school size £. The data obey a power law with 
exponent 7 = 0.16 ± 0.05. Using (Mj) and this value of 7, we 
obtain an independent estimate (3 = 0.25 ± 0.05. 



